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Abstract. Absorption imaging and spectroscopy can probe the dynamics of an 
ultracold neutral plasma during the first few microseconds after its creation. 
Quantitative analysis of the data, however, is complicated by the inhomogeneous 
density distribution, expansion of the plasma, and possible lack of global thermal 
equilibrium for the ions. In this article we describe methods for addressing these 
issues. Using simple assumptions about the underlying temperature distribution and 
ion motion, the Doppler-broadened absorption spectrum obtained from plasma images 
can be related to the average temperature in the plasma. 
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1. Introduction 

Ultracold neutral plasmas l a , formed by photoionizing laser-cooled atoms near the 
ionization threshold, stretch the boundaries of traditional neutral plasma physics. The 
electron temperature in these plasmas is from 1-1000K and the ion temperature is around 
1 K. The density can be as high as 10 10 cm -3 . Fundamental interest in these systems 
stems from the possibility of creating strongly-coupled plasmas 0, but collective modes 
recombination and thermalization [H| have also been studied. 

Charged particle detection techniques have traditionally been used for these 
experiments. However, optical absorption imaging and spectroscopy, demonstrated 
[5] using the Sr + 2 Si/ 2 — > 2 Pi/2 transition in a strontium plasma, opens many new 
possibilities. Images depict the density profile of the plasma, and the Doppler-broadened 
absorption spectrum measures the ion velocity distribution. Both can probe ion 
dynamics with 50 ns resolution. 

Qualitative interpretation of the images and spectrum is straightforward, but 
quantitative analysis is complicated by the inhomogeneous density distribution of the 
plasma, plasma expansion, and the possible lack of global thermal equilibrium for the 
ions. In order to address these effects, we are forced to make some simple assumptions 
about the form of the ion temperature distribution and expansion of the plasma. In 
this article we motivate these assumptions and describe our methods for quantitatively 
analyzing the data. 

Section El provides an overview of the plasma creation and absorption imaging 
technique. Section|3]describes the dynamics of the ions during the first few microseconds 
after photoionization. This provides the basis for understanding Sec. 01 which explains 
how the absorption spectrum is extracted from the images and how it relates to the 
temperature distribution. 

2. Experimental Overview 

The production of a strontium plasma starts with atoms that are cooled and confined 
in a magneto-optical trap (MOT) (Figs. QandEJ). This aspect of the experiment was 
described in jH]. The neutral atom cloud is characterized by a temperature of a few 
mK and a density distribution given by n(r) = noexp(— r 2 /2a 2 ), with a ~ 0.6 mm and 
n ~ 6x 10 10 cm -3 . The number of trapped atoms is typically 2 x 10 8 . These parameters 
can be adjusted. In particular, turning off the trap and allowing the cloud to expand 
yields larger samples with lower densities. 

To form the plasma, the MOT magnets are turned off and atoms are ionized with 
photons from the cooling laser and from a 10 ns pulsed dye laser whose wavelength 
is tuned just above the ionization continuum (Fig. [TJ). Up to 30% of the neutral 
atoms are ionized, producing plasmas with a peak electron and ion density as high 
as n 0e pa n 0i pa 2 x 10 10 cm -3 . The density profiles, n e (r) pa 7ij(r), follow the Gaussian 
shape of the neutral atom cloud. 
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Figure 1. Strontium atomic and ionic energy levels involved in the experiment, with 
decay rates. (A) Neutral atoms are laser cooled and trapped in a magneto-optical trap 
(MOT) operating on the 1 S'o — 1 P\ transition at 460.9 nm, as described in Atoms 
excited to the 1 P\ level by the MOT lasers are ionized by photons from a laser at 
~ 412 nm. (B) Ions are imaged using the 2 Si/2 — 2 P/2 transition at 421.7 nm. 2 P\/2 
ions decay to the 2 D 3 / 2 state 7% of the time, after which they cease to interact with 
the imaging beam. This does not complicate the experiment because ions typically 
scatter fewer than one photon during the time the imaging beam is on. 
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Figure 2. From [5]. Experimental schematic for strontium plasma experiments. The 
MOT for neutral atoms consists of a pair of anti-Hclmholtz magnetic coils and 6 laser- 
cooling beams. Atoms from a Zeeman-slowed atomic beam enter the MOT region and 
are trapped. 1 P\ atoms are ionized by the photoionizing laser. The imaging beam 
passes through the plasma and falls on a CCD camera. 



Because of the small electron-ion mass ratio, the electrons have an initial kinetic 
energy approximately equal to the difference between the photon energy and the 
ionization potential, typically between 1 and 1000 K. The initial kinetic energy for the 
ions is close to the kinetic energy of neutral atoms in the MOT. As we will discuss below, 
the resulting non-equilibrium plasma evolves rapidly. 

To record an absorption image of the plasma, a collimated laser beam, tuned near 
resonance with the principle transition in the ions, illuminates the plasma and falls on 
an image intensified CCD camera. The ions scatter photons out of the laser beam and 
create a shadow that is recorded by an intensified CCD camera. The optical depth 
{OD) is defined in terms of the image intensity without (hackground) and with (I p / asmo ) 



Absorption Imaging and Spectroscopy of Ultracold Neutral Plasmas 



4 




Y(mm) x(mm) 



Figure 3. Optical depth of an ultracold neutral plasma. The delay between the 
formation of the plasma and image exposure is 85 ns. The plasma contains 7 x 10 7 
ions and the initial peak ion density is n$i — 2 x 10 10 cm' 3 . Resolution is about 65 /^m, 
limited by pixel averaging to improve the signal-to-noise ratio. 

the plasma present, 

OD(x,y) = ln(I background (x,y)/I plasma (x,y)). (1) 

Figure El shows a typical absorption image. Section 0] describes how detailed information 
about the plasma is extracted from the optical depth. 

3. Ion Dynamics 

In order to understand the details of the image analysis, it is necessary to understand 
the dynamics of the plasma. The imaging probe is most sensitive to the ion dynamics, so 
we will concentrate on this topic. The behavior of electrons was studied experimentally 
in HIEIIII and theoretically in 00 El- 

Ions are created with very little kinetic energy, but their initial spatially 
uncorrelated state possesses significant Coulomb potential energy compared to the 
regular lattice that represents the ground state of the system [TQJIII]. As ions equilibrate 
and correlations develop, the kinetic energy increases. This process is called disorder- 
induced heating, or correlation heating, and it has been discussed in many theoretical 
papers. Early interest was generated by non-equilibrium plasmas created by fast-pulse 
laser irradiation of solid targets, for example O El EH] , and the problem has been 
revisited in the context of ultracold neutral plasmas |H1 HE1 HI] ■ 

Qualitatively, one expects the ion temperature after equilibration to be on the order 
of the Coulomb interaction energy between neighboring ions. A quantitative analysis 
|16j . assuming complete initial disorder and incorporating the screening effects of the 
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electrons, predicts an ion temperature of 
2 e 2 , ~ k , 




(2) 



3 47T£o&^B 2 



Here, k = a/Ac where Ad = (eokBT e /n e e 2 ) 1 / 2 is the Debye length. The quantity 
£7 = 7v e 2 /47r£ a * s ^h e excess potential energy per particle in units of e 2 /4:iTE a, where 
a = (47rnj/3)~ 1 / 3 is the Wigner-Seitz radius, or interparticle distance. iV, is the number 
of ions. U has been studied with molecular dynamics simulations ^Hj for a homogeneous 
system of particles interacting through a Yukawa potential, <p(r) = 4 ^ eor exp(— r/Ap), 
which describes ions in the background of weakly coupled electrons f. 

For typical strontium plasmas discussed here, k ~ 0.1 — 1, and Ad ~ 2 — 8/im. 
C7 ranges from —0.6 to —0.8, so Tj is close to e 2 /4ne akB as expected, k is related to 
the Coulomb coupling parameter for electrons, T e , through k = \/3T e . A system is 
strongly coupled when T > 1 [2.. r e « 0.1 — 0.5 for the systems studied here, so the 
electrons are not strongly coupled. This avoids excessive complications that arise when 
T e approaches or initially exceeds unity, such as screening of the ion interaction [7] , and 
rapid collisional recombination and heating of the electrons El El Ej , although we do 
see some signs of these effects, even in this regime. The ions typically equilibrate with 
Tj ~ 1 K, which gives Tj « 3, so the ions are strongly coupled. 

The time scale for disorder-induced heating is the inverse of the ion plasma 
oscillation frequency, l/uj pi = ^Jm^eo/n^e 2 , which is on the order of 100 nanoseconds. 
Physically, this is the time for an ion to move about an interparticle spacing when 
accelerated by a typical Coulomb force of e 2 /47T£oa 2 . This time scale is also evident in 
molecular dynamics simulations of ion-ion thermalization [7J |H1 El El El El El E] ■ 

It is interesting to note that under usual conditions in weakly interacting plasmas or 
even atomic gases, the time scale for relaxation of the two-particle distribution function, 
which describes correlations, is much faster than the collision time that governs the 
relaxation of the one-particle distribution function to the Maxwell-Boltzmann form. 
This is known as Bogoliubov's hypothesis [201 - For strongly-coupled plasmas, however, 
these time scales both become equal to the inverse of the plasma oscillation frequency 



As the two-particle distribution function equilibrates, the kinetic energy of the 
ions exhibits strongly damped oscillations at twice the ion plasma oscillation frequency. 
Intuitively, this can be understood as the oscillation of each ion in its local potential 
energy well. It is questionable whether this should be called an ion plasma oscillation or 
not because there is probably no collective or long range coherence to the motion. This 
behavior has been observed in molecular dynamics simulations of equilibrating strongly- 
coupled systems [T31 El El- The damping time for the oscillations is approximately 
n/ujpi for T > 5. Averaging over the entire density distribution, as we do in the analysis 
described here, obscures the oscillations. A different approach, which resolves regions of 

| As the number of electrons per Debye sphere (k~ 3 ) approaches unity, the Yukawa interaction ceases 
to accurately describe ion-ion interactions. For strontium plasmas studied here, this situation only 
occurs for the highest n e and lowest T e . It will be interesting to test Eq. 0for these conditions. 



B3|. 
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different density in the plasma, is required to clearly observe the oscillations, and this 
phenomenon will not be discussed further in this paper. 

For tdeiay > ^/upi, the ions have equilibrated and the thermal energy of the 
electrons begins to dominate the evolution of the plasma. Electrons contained in the 
potential created by the ions exert a pressure on the ions that causes an outward radial 
acceleration. This was studied experimentally in |3, and theoretically by a variety of 
means in [§]. The experiments measured the final velocity that the ions acquired, which 



was approximately v te rminai ~ yE e /mi. With the imaging probe, we now observe the 
expansion dynamics at much earlier times during the acceleration phase. 

As discussed in jH] and jS], a hydrodynamic model, which describes the plasma 
on length scales larger than Xd, shows that the pressure of the electron gas drives the 
expansion through an average force per ion of 

- -V(n e (r)k B T e ) ^rk B T ( 



- r (3) 

niir) af 

where the ion and electron density distributions are n e (r) « n»(r) = n 0i exp(— r 2 /2o~f). 
We assume thermal equilibrium for the electrons throughout the cloud 0. 
The force leads to an average radial expansion velocity for the ions, 

v(r,t de i ay ) = r rkBT * t delay . (4) 



The velocity is correlated with position and increases linearly with time. This does not 
represent an increase in the random thermal velocity spread or temperature of the ions. 
Due to the large mass difference, thermalization of ions and electrons [7j is slow and 
occurs on a millisecond time scale. 

Equation 0] for the average ion velocity assumes a constant electron temperature. 
Actually, as the plasma expands, electrons will cool. This can be thought of in 
terms of energy conservation or adiabatic expansion. It is possible to describe the 
expansion with a Vlasov equation that includes the changing electron temperature. 
For an initial Gaussian density distribution, the equations can be solved analytically 
and the expansion preserves the Gaussian shape with a 1/y/e density radius given by 
of(i) ~ of (0) + [kBT e (0)/mi]t 2 [HJI2I1I221- The experiments involving absorption imaging 
of the plasma, however, have concentrated on the first few microseconds of the expansion 
when the plasma size and electron temperature have not changed significantly. Thus we 
can safely use Eq. |3J The effects of the expansion are evident in the radial velocity that 
manifests itself in Doppler broadening of the ion absorption spectrum. 



4. Doppler-Broadened Spectrum 

To obtain quantitative information from the plasma images, we relate the OD (Eq.[IJ to 
underlying physical parameters. Following Beer's law, the OD for a laser propagating 
along the z axis is 

OD(x,y)= dzni(x,y,z)a[u,Ti(r)}, (5) 
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Figure 4. Absorption spectra of ultracold neutral plasmas. We plot the integral of 
the optical depth (Eq. ||JJ). The frequency is with respect to a Doppler-free absorption 
feature in a strontium discharge cell. Both spectra correspond to T e — 56 K and an 
initial peak plasma density of rtoi = 2 x 10 10 cm -3 . Data are fit with Voigt profiles, 
and the increase in linewidth for longer t^eiay is clear. 



where rii(x,y,z) is the ion density, and a[u, T»(r)] is the ion absorption cross section at 
the image beam frequency, v. The absorption cross section is a function of temperature 
due to Doppler broadening, and since we expect the temperature to vary with density, 
we allow a to vary with position. If we now integrate over x and y, or, in reality, sum 
over the image pixels multiplied by the pixel area, we get the spectrum 

S{y) = J dxdyOD(x,y) = J d 3 r rii(r)a[u,Ti(r)], (6) 

as a function of the image laser detuning §. As we vary the detuning, we obtain 
absorption spectra as shown in Fig. |3J The rest of the paper will deal with the 
relationship between spectra such as these and the underlying temperature distributions 
of the ions. 

The absorption cross section for ions in a region described by a temperature T«, is 
given by the Voigt profile 

a(u T) = I ds 3 * nX2 - - e -(s-"o) 2 /2« D ( Tl r m 



where aniTi) = ykBTi/rrii/X is the Doppler width, and 7 e // = "Jo + liaser is the effective 
Lorentizian linewidth due to the natural linewidth of the transition, 70 = 2n x 22 x 10 6 
rad/s, and the laser linewidth, 7; aser = 2n x (10 ± 2) x 10 6 rad/s. The center frequency 
of the transition is uq = c/A, where A = 422 nm. The "three-star" symbol, 3* = 1, 

§ We can also fit OD(x,y) to a two dimensional Gaussian, as described in 0, and identify 
J dxdyOD(x,y) w 2Tr<7i X ai y OD max , where (T,- x and <7j a are the transverse sizes of the absorption 
profile, and OD max is the peak optical depth. This sometimes has signal-to-noise ratio advantages 
over integrating the entire image, but both approaches should give the same result. 
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accounts for the equal distribution of ions in the doubly degenerate ground state and 
the linear polarization of the imaging light [23 . 

Extracting the ion temperature from the spectrum is complicated by the fact 
that we do not expect to have global thermal equilibrium of the ions during the 
first microsecond after plasma formation. Global thermal equilibrium will occur on 
a hydrodynamic time scale, Oijv , which is on the order of ten microseconds, where v is 
the ion acoustic wave velocity ||. 

We need a method to relate the Doppler broadening of the spectrum to a 
characteristic temperature. We will show that for reasonable models of the ion 
temperature distribution, the temperature extracted from a fit of a Voigt profile to 
the spectrum yields a temperature that is a close approximation of the average ion 
temperature in the plasma. More thorough checks of this assignment, such as with 
molecular dynamics simulations, would be valuable, but this serves as a useful working 
definition for discussing the data. 

For times longer than n/upi, but short compared to hydrodynamic times, we expect 
local thermal equilibrium at a temperature approximately given by Eq. |21 Under this 
model, the temperature varies across the plasma because the interparticle spacing, a, 
varies with density. Given T e and n^, it is possible to use Eq. |2]in an iterative recipe 
to find the local ion temperature. This procedure utilizes an expression [21] for the 
excess potential energy U in terms of and k that is valid for k < 5. This yields ion 
temperature distributions as shown in Fig. 

To simplify the discussion, we first neglect electron screening. This accurately 
describes the plasma for k<1 (high T e and low n e ). To a good approximation, we can 
then say 

Ur) = C C ] = T hmax e- r2 ^, (8) 

where T i)max is the temperature corresponding to the peak density at r = 0. Molecular 
dynamics simulations of ultracold plasmas that neglect electron screening display 
this thermal distribution when the ions have equilibrated locally (t de i ay > n/u) pi ) [25J. 
C = 0.45 is a constant determined by numerical evaluation of the parameters in Eq. El 
for small k \ The average ion temperature, assuming Eq. |HJ is 



Combining the expressions for the integral over the optical depth (Eq. IHJ), and the 
expression for the absorption cross section (Eq. |7J) yields 

S(u)= [d 3 r ni (r) [ds^^ } - N9 _ 1 e -(s-u )V2a D m 

J J 2 1 + 4(^P v& D [7Kr)] 

|| In principle, it is possible to obtain spectra from small regions of the cloud, so as to examine the 
local dynamics and avoid averaging over the plasma. Developing this capability will be the subject of 
future studies. 

f As k increases, Eq. 03 still describes the ion temperature if C decreases and becomes a function of 
position in the plasma, as shown in Fig. El 
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Figure 5. Ion temperature profiles. Peak densities (no,i ~ n o,e) and electron 
temperatures are given in the legend. Screening refers to the iterative calculation 
of temperature using Eq. [21 No screening refers to Eq. HO which neglects electron 
screening, and is thus independent of electron temperature. In the limit of low density 
and high electron temperature, both expressions give the same result. 



If we insert (Eq. |SJ) for Tj(r), we find 



-e^exp[-e^( S - ^ )72< m J. (11) 



Here aD )max is the Doppler width corresponding to T i>max . As shown in Fig. El the 
average over the density distribution of the Doppler factor can be replaced by a single 
Doppler distribution with temperature T = 0.59T imax = 0.91T ijave . 



27TCT Di „ 

^— exp[-(s-u ) 2 /2d 2 D }, (12) 



where cr^ = ^0.91kBT itave /mi/X. So if the velocity distribution of the ions is 
characterized by local thermal equilibrium of the form in Eq. |H1 fitting the spectrum, 
S(u), to a Voigt profile yields a good approximate measure of the quantity Ti >ave . 

If we want to include the effects of electron screening, we can follow exactly the 
same procedure, except in our model of the spectrum (Eq. HQJ) we use the iterative 
recipe for ion temperature using Eq. El instead of the simple Eq. |H1 This slightly 
more complicated numerical analysis shows that if Eq. El describes the temperature 
distribution, then the temperature extracted from the Voigt profile is betwen 0.91 and 
0.95 of Ti >ave . In this case, the average ion temperature is found by averaging Eq. |2]over 
the density distribution. 
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Figure 6. Assuming Ti(r) is given by Eq. |SJ the Doppler profile averaged over 
the density distribution can be replaced by a Doppler profile characterized by a 
single uniform temperature, T = 0.59Ti_ max = 0.91Ti jOt , e . The statistical weight 
is proportional to the probability of finding ions with that Doppler shift, and the 
horizontal axis is the Doppler shift in units of &D,max' 



If we have global thermal equilibrium, which we expect for long t delay, then the 
temperature is uniform, and, by definition, the Voigt profile will yield the average 
temperature. 

The expansion of the plasma, Eq. HI also contributes to the velocity distribution. So 
we must expand our model. Fortunately, the expansion is well understood [SI El 122] and 
determined by the parameters T e and cij. Independent analysis of the images determines 
(7j. T e is approximately equal to 2E e /3kB, where E e is the detuning of the photoionizing 
laser above resonance [T], although for high density and low electron temperature, we 
can expect a small increase of T e above 2E e /3kB due to electron heating from three-body 
recombination [9 , disorder-induced heating [Zj, and continuum lowering jH]. 

We will first treat the expansion using the simple model of ion temperature, Eq. 
IH1 The expansion velocity profile gives rise to an average Doppler shift of the resonant 
frequency that varies with position, 

r- / x TksTet delay n 
OV = V z A = » cosU 

miafX 

= 0- D ,ave—-^-COs9, (13) 

where rcosO is the displacement from the center of the cloud along the direction of laser 
propagation. We have introduced a characteristic time for the expansion, 




*exp = ^/r-^ry-^-' ( 14 ) 

which is the time at which the Doppler shift due to expansion, at rcosO = o~i, equals 
o~D,ave- In other words, this is the time at which the Doppler broadening due to expansion 
becomes comparable to the thermal Doppler broadening. For typical plasma conditions, 
such as T e = 30 K, T itave = 1 K, and (Tj = 1 mm, Eq. H3J yields t exp = 3.4 /xs. 
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Figure 7. Doppler profile averaged over the density distribution as in Fig. El 
but now including the effects of expansion. The dashed or dash-dot numerical 
expressions are full integrals over volume shown in Eq. 1151 The solid curves are the 
Gaussians for Doppler broadening arising from a single temperature T = 0.95 Ti^ ave [1 + 
(1/0.95) (t/t exp ) 2 }, where t exp is the characteristic time scale for the expansion described 
in the text. At all times, the single Doppler profile is a good approximation to the full 
integral. 



We include the position-dependent Doppler shift in the exponent of the Gaussian 
describing the Doppler broadening in the Voigt convolution (Eq. fTT| . which becomes 



/ 



j3 -r 2 /2a 2 ^ -r 2 /12cr 2 r r 2 /6<r 2 / T ^delay n\2 /n 2 1 

d rn 0i e 1 » -= e 1 »exp[-e 1 > (s - u - a D ,ave — — —cosO) /2<J D ma J 

V 2,TlO~£i max 0~i t exp 

exp[- n 2 , m d- (15) 



V2no- D (T i>eff ) 2a%(T iieff) 

In the second line of Eq. [Hi we have given the Doppler profile, corresponding to a 
temperature of 

T, e// = CT^ e [l + i(^) 2 ], (16) 

that approximately equals the original full integral, including variation in temperature 
across the plasma and effects of expansion. Here, C = (0.95 ± 0.05). Figure shows 
that the agreement is good over the entire time scale of the expansion. 

If we include screening effects and use Eq. Elto find the temperature, then numerical 
simulation shows that the extracted temperature is also given by Eq. EE with C ranging 
from 0.91 to 0.95. If the plasma has reached global thermal equilibrium, then it can be 
shown analytically that T = T ijave [l + (t de i ay /t exp ) 2 ] exactly. This suggests that when 
fitting a Voigt profile to the spectrum derived from the integral over the optical depth, 



2 l + 4(^) ! ^„(T,, t „, 

the extracted effective temperature can be related to the average temperature in the 
plasma through Eq. One can show analytically that Eq. ^] has the correct form 
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Figure 8. Evolution of the ion temperature in an ultracold neutral plasma. The initial 
peak ion density is noi = 5 x 10 9 cm -3 and Oj = 780 /im. The plasma contains 6 x 10 7 
ions. (A) The effective ion temperature corresponding to the Doppler broadening 
extracted from Voigt profile fits of data such as in Fig. 0J The solid line (Eq. I16(l is the 
effective temperature one would expect for Ti^ ave — 0.9 K and T e = 2£ , e /3fcs = 40 K. 
Equation^|includes the effects of expansion but neglects early equilibration. (B) The 
average ion temperature is obtained by dividing by 0.95 [f + (l/0.95)(td e iay/t e xp) 2 ]- 

in the limit when all broadening is negligible except the Doppler broadening due to 
expansion. 

Our motivation of Eq. El relies on certain models of the ion temperature. They are 
reasonable approximations to the distributions seen in molecular dynamics simulations 
except at short times, t < n/uj pi « 500 ns, during the initial rapid disorder-induced 
heating and oscillatory phase when the kinetic energy distribution as a function of 
temperature is not well known. Because this time scale varies with density, we 
expect different regions to equilibrate and oscillate at different times. Extracted ion 
temperatures should be viewed as more qualitative for these early times. 

Equation is the central result of this paper. It enables us to analyze the data 
in a clear and well-defined way, as shown in Fig. |Hl It also separates the various 
contributions to the Doppler broadening of the spectrum. For example, Fig. EK shows 
the time dependence of T^ e ff for a typical plasma. Figure EB shows the extracted 
Ti t ave with the expansion effects removed. The only fit parameter in the model line 
in Fig. ER is T i (M , e . This procedure allows us to determine T iAve to an accuracy of 
about 0.2 K. Alternatively, after removing the effect of expansion by scaling the data by 
0.95 [1 + (1 / 0. 95) (t delay /texp) 2 ], as in IH)3, the equilibration temperature can be identified 
by the plateau. Scaling the data also improves the visibility of the disorder-induced 
heating during the first l/u> pi w 180 ns. When there is significant electron heating 
[3 HE IE], T e will exceed 3E e /2kB, and T e then becomes a fit parameter that provides 
information on the electron dynamics. 

Using this recipe for analyzing the data we can conduct quantitative studies of the 
rate and equilibration temperature for disorder-induced heating. We can also examine 
the factors affecting ion temperature oscillations, electron temperature, and plasma 
expansion. These will be the subjects of future work. 
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